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Abstract 



The physics of a parity anomaly, potentially observable in a narrow-gap semi- 
conductor, is revisited. Fradkin, Dagotto, and Boyanovsky have suggested 
that a Hall current of anomalous parity can be induced by a Peierls distortion 
on a domain wall. I argue that a perturbation inducing the parity anomaly 
must break the time reversal symmetry, which rules out the Peierls distortion 
as a potential cause. I list all possible perturbations that can generate the 
anomaly. 

I. INTRODUCTION 

Condensed-matter and particle physics are cross-fertilizing fields. For instance, exclusion 
of magnetic flux from a superconductor was interpreted by Anderson as generation of mass 
for a gauge boson.El Popularized in particle physics by Higgsp this effect plays a central role 
in the electroweak theory. In the other direction, statistical (Chern-Simons) gauge fields, 
a field-theorist's toy, facilitate better understanding of the quantum Hall effect.u A closely 
related phenomenon, parity anomaly in (2 + l)-dimensional electrodynamics, still awaits its 
experimental discovery in condensed matter. 

In a nutshell, the effect is rather simple. Two-component Dirac fermions with a mass m 
react in a peculiar way to an external magnetic field B. The symmetry between occupied 
(Fermi-sea) and empty fermion states is violated and a nonzero fermion density appears in 
the vacuum state. The density of fermions — each carrying charge ±e/2 — is such that the 
average flux per particle is hc/e: 

_Q_ m_^B^ , . 

^~ A~ \m\2hc' ^ ^ 

By relativistic invariance, application of an electric field in the plane induces a Hall current 
with quantized conductivity: 
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J^ = {(^/c)eijEj, a = -r^^- (2) 



Both the induction of charge by magnetic field (P and the appearance of a Chern-Simons 
current (^ violate the symmetries of parity and time reversal, hence the name "parity 
anomaly". In certain materials, electrons behave as Dirac fermions, albeit with a small 
"speed of light" , so it seems natural to look for the anomaly in condensed matter. 

On a deeper level, there is a subtle problem, known as fermion doubling, which often 
prevents the anomaly. For example, electron states in a sheet of graphite are well described 
at low energies as 2-component massless Dirac fermions in 2+1 dimensions. Because the 
unit cell contains 2 lattice sites, there are 2 fermion species living at inequivalent Fermi 
points in the Brillouin zone. Semenoff observed that a symmetry breaking mass term m will 
be induced if the inequivalent sites are populated by different atoms. Unfortunately, the 2 
species have mass terms of opposite signs.a The total anomalous current @ vanishes. 

Later, Fradkin, Dagotto and BoyanovskyO (FDB) suggested a possible way to circumvent 
the problem of fermion doubling. Lead chalcogenides PbTe, PbSe, and PbS are semiconduc- 
tors with a narrow gap (0.15-0.3 eV) between conduction and valence bands.Q Low-energy 
fermion quasiparticles are concentrated around 4 inequivalent L points in the Brillouin zone, 
±(7r/a)(l, 1, 1) and three others related by cubic symmetry. The band structure near each 
L point is such that the quasiparticles resemble 4-component massive Dirac fermions (in 
3+1 dimensions). 4 components result from 2 inequivalent sites in a unit cell and 2 pos- 
sible projections of a spin 1/2. FDB pointed out that a stack fault in a crystal creates a 
domain wall, on which one finds 2+1 dimensional massless Dirac fermions. In this setting, 
certain perturbations can induce a mass term of the same sign in all 4 fermion species. The 
anomalous current (Q) does not have to vanish. 

In a more detailed paper ,0 FDB suggested that a symmetry-breaking mass term m can be 
induced by a Peierls distortion, which is present in some lead chalcogenides.Q That does not 
sound right: a Peierls distortion violates parity but leaves time reversal intact, therefore it 
cannot possibly induce a Chern-Simons current @. Haldane has argued that this incarnation 
of parity anomaly is caused by an uuDhysical lattice Hamiltonian used by FDB, which is (or 
seems to be) odd under time reversaLEl Upon further refiection, this argument does not work: 
if a continuum limit derives from a model with lattice potential and spin-orbit interaction, 
it cannot violate the time-reversal symmetry. Then why does the FDB Hamiltonian couple 
spin and linear momentum via the term a-k — that seems to break T? The answer: the a 
matrices here are not the electron spin variables. Symmetry under time reversal need not 
be broken. The existence of the anomalous Chern-Simons current is then questionable. 

Being field theorists, FDB didn't bother to derive their Hamiltonian (nor its continuum 
version) from any physical model of electrons in a lattice potential and with the spin- 
orbit interaction. Their phenomenological Hamiltonian was merely tailored to correctly 
reproduce the energy spectrum of the conduction and valence bands near the Fermi level. 
Therefore the a matrices in it have nothing to do with the actual electron spin. But without 
knowing which physical variables the a matrices represent, one cannot learn how this or 
that physical perturbation (e.g., a Peierls distortion) couples to the fermion zero modes. 
Then it is impossible to determine correctly the sign of charge — or the direction of Hall 
current — induced by a symmetry-breaking perturbation. 

To clear up the matter, I have derived a correct continuum approximation for low-lying 
electron states both by appealing to symmetry arguments (the easy way) and by using the 
standard k-p method combined with a tight-binding approximationciEJ (see Appendix). I 



have studied the behavior of massless fermions bound to a domain wall and enumerated all 
symmetry-breaking perturbations that raise or lower the zero modes in magnetic field, thus 
inducing surface charge (^. The results of this study are unambiguous. Already from Eq. 
m it is clear that m should be odd under both parity and time reversal (a pseudoscalar) . 
A Peierls distortion therefore will not work because m = q-n is a genuine scalar (q is the 
vector of atomic displacement and fi is the normal to the domain wall). On the other hand, 
certain kinds of magnetic order may induce an anomalous Hall current. 

II. LOW-ENERGY HAMILTONIAN 

The FDB Hamiltonian (in the continuum limit) can be derived from a tight-binding 
model of noninteracting electrons that includes electron kinetic energy, periodic lattice po- 
tential, and spin-orbit interaction:ll2l 

This Hamiltonian is invariant under parity and time reversal. 

The triplet of Pauli matrices a denotes the spin operators. Introduce operators of isospin 
Ti, T2, T3 to describe the two sublattices in the rocksalt structure. Eigenvalues T3 = ±1 corre- 
spond to lead and chalcogen sites. Near each of the four L points, lowest-energy states have a 
4-component wavefunction (2 spin components x 2 sublattices) and resemble Dirac fermions. 
By symmetry arguments alone, one can anticipate the correct form of the Hamiltonian in 
the continuum limit. 

Define the operations of parity and time reversal as 

Pip{r) = -^(-r), Tip{r) = cTyi;*{r). (4) 

The minus sign accounts for odd parity of p orbitals. Near an L point in the Brillouin zone, 
e.g., p = (7r/a)(l, 1, 1), the one-particle Hamiltonian of a Dirac fermion is limited to the 
following P and T-invariant terms: 

H = -ivTi{p ■ V) - w\T2{[a X p] • V) + t^Mv'^. (5) 

Here p = p/p is odd under both parity and time reversal (directions p and — p describe the 
same L point); A is the relative strength of the spin-orbit interaction, and f is a fermion 
velocity. For simplicity, I set A = 1, which gives a spherically symmetric energy spectrum 
near the L points. The same form of H is obtained from a tight-binding model, see Appendix 
^ The low-energy Hamiltonian thus has a Dirac form H = —ia-'V + jSM. The standard 
Dirac matrices have the following representation: 

a = Tip + T2[a X p], ^ = 70 = r3, 

-ij = -iPd = T2P - Ti[a X p], 'j^ = ~Ti{a-p), (6) 

E = -7^a = ((T-p)p + Tspx [(Txp], ■i7°7^ = T2{(T-p) 



III. MASSLESS FERMIONS ON A DOMAIN WALL 

Consider a stack fault (Pb ^^ Te) in a plane perpendicular to a unit vector h. Then M 
is a function of X|| = r-fi. More specifically, the normal to the domain wall n points in the 
direction of increasing M: 

M{xii) -^ ±|Mo| for a;|| = r-n -^ ±00, (7) 

Energy eigenstates bound to the domain wall can be written in the form '?/'(r) = u{x\\)ip{'x±), 
where x^ = nx[rxn] are coordinates within the plane. The scalar u{x\\) and bispinor V'(x_|_) 
satisfy the following equations:0 

u'{xii) = ~M{xii)u{xii), (8) 

(-^7•n)^(x^)=V'(xx), (9) 

-wa-V±ip{x±) = E^{^i_) (10) 

The first two equations (H-^ indicate that these states are fermion zero modes in the 
direction perpendicular to the wall.li3 The last two (PHT0|) describe 2-component massless 
fermions in 2+1 dimensions with the Hamiltonian 

H = vd-{-iV^--A^). (11) 

Application of magnetic field B = Bri perpendicular to the domain wall results in a 
spectrum with Landau levels:t3 

H^ = v^i-iV^--AA S-n. (12) 

The energy spectrum is symmetric with respect to charge conjugation: 

E = 0, ±^;(2|eS|/c)^/^ ±t;(4|e5|/c)^/' . . . (13) 

Zero modes {E = 0) are eigenstates of the "spin" component S ■ n with the eigenvalue 
sgn(ei?). All Landau levels have orbital degeneracy \eB\/2TTc per unit area. 

In a half- filled system (no dopants), the zero modes are exactly at the Fermi level and 
thus have occupation numbers 1/2. The energy spectrum is particle-hole symmetric and 
the domain wall is not charged. An arbitrarily small perturbation can shift the zero modes 
above or below the chemical potential (assuming it stays pinned at 0). With the particle- 
hole symmetry broken, the domain wall gets charged: each zero mode contributes charge 
+e/2 (empty) or — e/2 (filled) to the domain wall (e < 0).ll3 

Although a condensed-matter system lacks the true relativistic invariance and quantiza- 
tion of Hall conductivity (0) does not follow automatically from Eq. [l|, there is a thermody- 
namic identity that establishes this relation:li^ 



dp 
" = -95 



m e^ 

\m\ 2/1' 



(14) 



with the derivative taken at a constant chemical potential. This result is valid in the absence 
of low-energy excitations. 



IV. SYMMETRY-BREAKING PERTURBATIONS 

Imagine now that a uniform perturbation V (such as a Peierls distortion or the actual 
Zeeman term) is apphed to fermions on the domain walL In the presence of the Landau 
gap, the shift of the zero modes can be computed to first order in V using the standard 
perturbation theory. Because there is no spin or isospin degeneracy and V is uniform, 
AE = f^p^Vip d^r. 

Recall that ip is an eigenstate of 3 commuting matrices: 

(S-n)V^ = sgn(e5)^, (15) 

(-n-n)V = ^, (16) 

{ij^j^) ^ = sgn(eS) ^p. (17) 

(The third matrix 7*^7^ is simply the product of the first two.) To compute tp'^Vip, one can 
write V SiS a. superposition of 15 traceless Hermitian 4x4 matrices 

{H,} = {7O, 7^ ^7^^ «, -^7, S, 70s} : 

15 -, 

^^^^ = E jT^^VH,) ^%^ (18) 

Because operators Hi either commute or anticommute with one another, only 3 of them 



(p!5|-[T7D need to be included in the sum: if some Hj anticommutes with one of them, its 
expectation value in the state ip vanishes. 

Thus operators (pr^-[T7|) exhaust all the handles through which external perturbations 
can tickle the zero modes and possibly induce charge on a domain wall. Such a perturbation 
will contain F^-S, —iF^- 7, or 2^57^7^. The symmetry breaking term F^ is a pseudovector, 
F^ is a vector, and F^ is a pseudoscalar. The area density of charge induced by them on a 
domain wall (for a single fermion species) is 

,^ „, e^(B-n) ,^ ^.e\eB\ ,^,e^(B-n) ,,^, 

-^«""'-"'4rf^- -^^■'"'^■"'idr- -^«"'^^)4rf^- (''" 

It is now evident that a vector perturbation F^, such as a Peierls distortion, cannot 
induce a Chern-Simons current: the sign of induced charge is not sensitive to the direction 
of magnetic field. By the thermodynamic identity (|14D, there will be a normal Hall effect 
with current reversing the direction when magnetic field does. 

One obviously needs a pseudoscalar perturbation (odd under both P and T) in order to 
relate a scalar (charge) to the pseudoscalar B-fi. This is why a pseudoscalar F5 will work 
(and a vector F^ will not). Alternatively, one can use a pseudovector Fs (e.g., staggered 
field of an antiferromagnetic order parameter) to construct a pseudoscalar (Fs-n). 

To corroborate these general considerations, I discuss in some detail several symmetry- 
breaking terms that might exist in a physical system. 



A. Peierls distortion 

We need to add to the 3 + 1 dimensional Dirac Haniiltonian (j^) a term that breaks P 
but not T and does not affect electron spin. Such a term is 

V^ = r2(q-p) (20) 

(both T2 and p are T-odd). The vector q characterizes the direction and length of the Peierls 
distortion. It couples to the operator — 27-11 (plGl). With the aid of Eq. |1^, 

AE = ^Tr (r2(q-p)(-27-A)) = (q-p)(p-n). (21) 

The charge density is unchanged when magnetic field is reversed B ^ — B. By extension 
(0), the Hall conductivity is sensitive to the direction of B. For a wall perpendicular to q, 
summation over the four L points gives 

a = ^sgn(e5). (22) 

This describes a perfectly normal Hall current, not a Chern-Simons current (|^). No parity 
anomaly here. 



B. Zeeman effect: uniform magnetic field 

Consider Zeeman interaction V = —ge(B-cx)/2mc induced by uniform magnetic field B 
perpendicular to the domain wall. 



AE = -^ -Tr ((B.a)(S.n)) sgn(ei^) = -^(p-A)^- (23) 

At a constant chemical potential, 

Q 4e\eB\ 



A Anc 
Again, a normal Hall effect results: 



sgn(^). (24) 



a = -—sgn{geB). (25) 



C. Zeeman eff"ect: staggered magnetic field 

In a system with antiferromagnetic order, one expects staggered magnetic field Bgt to 
induce a Zeeman term V = —T3ge(Bsf(^)/2'mc. It couples to the same operator of Dirac 
"spin" (|15D: 



Ai? = -£^^Tr(r3(B,fa)(S.n))sgn(e5) = -^^|^^([px[B,tXp]].n). (26) 



In the simple case with Bgt = -Bstn, 

Q _ Ae^B 
A Aire 
and 



sgn{geB,t) (27) 



sgn(5). (30) 



a = -—sgn{geB,,). (28) 

The Hall current depends on the direction of the staggered field Bgt and not of the uniform 
field B. One can call this a parity anomaly, but clearly there is no magic involved: this is a 
Hall effect caused by the staggered magnetic field. 

D. A pseudoscalar? 

Finally, one can imagine coupling directly to the third conserved quantity (0). The only 
perturbation that will do is 

V = 5T2ia-p). (29) 

The sign of induced charge depends on the direction of magnetic field: 

Q _ Ae'^B 
A Aire 

Hall conductivity is insensitive to the sign of B: 

a = ^sgn{5). (31) 

This term definitely generates a Chern- Simons current. I must admit though that I do 
not understand what perturbation could cause it. 

V. CONCLUDING REMARKS 

I have presented a critique of the parity anomaly (Eqs. |l| and ^ suggested for 2+1 dimen- 
sional fermions in a narrow-gap semiconductor .B Both symmetry considerations and a correct 
continuum treatment of the model indicate that a parity anomaly requires a pseudoscalar 
[P and T-odd) or pseudovector (P-even and T-odd) symmetry-breaking term. Therefore, a 
Peierls distortion, which is a vector (P-odd and T-even), cannot generate a parity anomaly, 
contrary to the suggestion of FDB. 

Antiferromagnetic order, which violates the time-reversal symmetry, could lead to an 
observable Chern-Simons current in the system. Unfortunately, the prospects of observing 
such an effect are rather dim: not only would this require antiferromagnetism, but also the 
Neel order parameter must be (anti)parallel to the applied magnetic field. In a Heisenberg 
antiferromagnet, Bgt stays orthogonal to B, in which case the entire effect vanishes. 

In hindsight, this foray into the dreamworld of field theory is a reminder to a condensed- 
matter physicist: trust your intuition. Anomalous Hall current of Dirac fermions, whether 
directly in d = 2 + 1 dimensions or on a domain wall in rf = 3 + 1, is always caused by 
magnetic order present in the system.QBt3 
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APPENDIX A: DIRAC FERMIONS FROM A TIGHT-BINDING MODEL 

Lead and a chalcogen have 2 and 4 electrons in their 3 external p orbitals {s orbitals give 
deeply lying bands, which can be ignored). The p bands are thus half filled. The one-particle 
Hamiltonian including lattice potential and spin-orbit interaction is 

Omitting spin for clarity, one can write an energy eigenfunction in the Bloch form 

Px,Py,Pz Pb,Te 

^p(R) = J2 E i^'^'i'R - r - Ar'') (l)^\r)e'P-^'+^''"\ (A2) 

i a 

Here r labels unit cells in the rocksalt lattice, Ar° are the coordinates of site a within a 
unit cell and '?/'"*(R) is the wavefunction of the i-th orbital on sublattice a. We are to find 
the Bloch coefficients 0p*(r). With spin included, 0"*°^ is a 12-component wave function. 

It follows from symmetry considerations alonell2l that the two sublattices (lead and chalco- 
gen sites) are decoupled at the L points, such as p = ±(7r/a)(l, 1, 1). In this case, the 
Hamiltonian takes on a simple formO 

Ha = 2Wa cos (oL-p) + A„a-L, (A3) 

where the first term comes from the overlap of different p orbitals (between second neighbors) 
and the second, from the spin-orbit interaction. The eigenstates of the Hamiltonian ( |A3[ ) 
are Kramers doublets with J-p = ±1/2 or ±3/2. According to Volkov et a/.JliJthe spin-orbit 
coupling is the weakest perturbation and the two bands around the Fermi level are derived 
from states with L-p = 0. In this approximation, El the eigenstates are 




^--{r) = ^\ 1 |®|a, = a), (A4) 

in the basis oi Px,Py,Pz orbitals. 

To calculate the matrix elements of the Hamiltonian near an L point, at lattice momen- 
tum p ± k, I use the standard "k-p" method." Instead of shifting momentum away from a 
symmetry point p ^ p ± k in the wavefunctions ([A^), one transforms the Hamiltonian 



H ^ e-^'^-^i/e*^-^ = H+^k-[-zV+ ^""JJf^^ ) ■ (A5) 




The matrix elements of the second term "(k-p)/r7i" are then evaluated in the Hilbert space 
of the two Kramers doublets closest to the Fermi level at the L point. This is a good 
approximation for small k, provided that other states lie at energies much larger than the 
gap A = (Epb - ^Te)/2. 

In the tight-binding limit, overlap oi Px orbitals on different sublattices occurs along the 
X direction only. In the case oi px orbitals at momentum px = Tc/a (Fig. |I]), operators — zV 
and V\^(R) have the following nonzero matrix elements: 

{Fh,x = 0\{-idx)\Te,x = a/2) = Ci, (A6) 

(Pb, X = 0\dxV\Te, x = a/2) = 1C2, (A7) 

where Ci and C2 are some real constants. In the basis of states px, Py, Pz, the resulting k-p 
term reads 

/ A;,. \ / kyCTz - K(jy \ 

CiTi A;^ + CsTa A;,a, - A;,a, (A8) 

\OOA;^y \ kxCJy - kyax J 

near the L point p = (7r/a)(l, 1, 1). The operators of isospin r act according to 

nlPb) = |Te), r2|Pb) = z|Te), rglPb) = |Pb), TslTe) = -|Te). (A9) 

Matrix elements of the complete Hamiltonian between the states (^^) are: 

ICiTiikx + ky + k^) + \C2T2 [{ky - k^)(Tx + {kz - kx)(Ty + {kx - ky)a^] 

+ i(Epb-^Tc)r3. (AlO) 

This yields the anticipated continuum result (^ 

H = t;(k-p) Ti + Xv(k- [px a]) ra + Mv'^ T3, (All) 

where p is the unit vector in the direction of the L point p = (7r/a)(l, 1, 1). An identical 
result is obtained for the 3 other L points, which differ by the direction p. 

On a final note, the unit vector p = p/p is odd under parity and time reversal defined 
in Eq. ^. It can be made even by redefining them as 

P^ir) = -r^^Pi-r), T^{v) = r^cXyriv). (A12) 

This is the familiar standard representation of parity and time reversal in relativistic field 
theory.O 
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Te Pb Te 

FIG. 1. Calculation of overlap integrals for adjacent Px orbitals. Momentum px = vr/a. 
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